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ABSTRACT 

According to our common understandings, the original surface of a short- 
period comet nucleus has been lost by sublimation processes during its close 
approaches to the Sun. Sublimation results in the formation of a dust mantle on 
the retreated surface and in chemical differentiation of ices over tens or hundreds 
of meters below the mantle. In the course of NASA's Deep Impact mission, 
optical and infrared imaging observations of the ejecta plume were conducted by 
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several researchers, but their interpretations of the data came as a big surprise: 
(1) The nucleus of comet 9P/Tempel 1 is free of a dust mantle, but maintains its 
pristine crust of submicron-sized carbonaceous grains; (2) Primordial materials 
are accessible already at a depth of several tens of cm with abundant silicate 
grains of submicrometer sizes. In this study, we demonstrate that a standard 
model of cometary nuclei explains well available observational data: (1) A dust 
mantle with a thickness of ~ 1-2 m builds up on the surface, where compact 
aggregates larger than tens of micrometers dominate; (2) Large fluffy aggregates 
are embedded in chemically differentiated layers as well as in the deepest part of 
the nucleus with primordial materials. We conclude that the Deep Impact results 
do not need any peculiar view of a comet nucleus. 

Subject headings: comets: general — comets: individual(9P/Tempel 1) 



1. INTRODUCTION 



Short-period comets are originally composed of ices and dust that were present in the 
cold comet-formation region of the primordial solar nebula. When the comets leave the 
trans- Neptunian region and approach the Sun, volatile components of its nucleus begin to 
sublime. Expanding volatile gases drag light dust having a large surface-to-volume ratio and 
leave preferentially heavy dust of a small surface-to-vo lume ratio on the surface of a comet 
to form a dust mantle 



isee, e.g.. 



Prialnik et al.l |2004| . for a review). Indeed, optical and 



infrared observations of cometary comae suggest that large compact aggregates of tens of 
micrometers or larger lie in a dust mantle with depleted volatile materials, while large fluffy 



aggre gates are embedded with volatile materials below the dust mantle (IKolokolova et al. 



20071 ). Dust mantle formation itself does not differentiate the composition of the aggregates. 



while chemical differentiation of a nucleus takes place owing to preferential sublimation losses 
of highly volatile materials in deeper parts of the nucleus. Pristine amorphous water ices 
also crystallize to form a crystalline ice crust just below the dust mantle. In additi on, impact 
cratering by interplanetary bodies affects the surface evolution of the nucleus (e.g.. [Fernandez 



198ll ). Therefore, primordial materials in comet nuclei remain only in the deepest region far 
below the dust mantle. The above mentioned scenario for the evolution of a comet nucleus 
(hereafter referred to as the standard model) is applicable to comet 9P/Tempel 1 (hereafter, 
Tl), the target of NASA's Deep Impact (DI) mission: Numerical studies on thermal evolution 
of the Tl nucleus predict the formation of a dust mantle with a thic kness of ~ 1 m and a 



crystalline ice crust w ith a thickness of 40-240 m below the mantle (IBar-Nun et al.l Il989 



De Sanctis et al.l 120071 ). 
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Observations of Tl during the DI event gave us a good opportunity of examining the 
standard m odel for the Tl nucle us, because DI excavated the Tl nucleus up to ~ 10-20 m 
depth (e.g.. ISunshine et al.ll2007l ). Plenty of DI observational results support the standard 



model: (1) The presence of a dust mantle on T l is evidenced by a low thermal inertia 
of the nucleus estimated from its thermal maps jA'Hearn et allboosh : (2) Near-infrared 
spectroscopic observations of volatil es excavated by the D I event revealed thermal processing 
of the layers below the dust mantle ( Mumma et al.|[2005l ): (3) The existence of large compact 
dust in the pre-impact Tl coma has been re ported by many researchers (jSykes fc Walker 
19921 : iFarnham et al.l 120071 : iReach et al.l 120071 ) . The existence of such large dust in the pre- 
impact Tl is consistent with the standard model, because the dust mantle formation depletes 
only fluffy aggregates. 

However, recent several researchers drew a peculiar picture of a comet nucleus from 
their observations of the ejecta plume pr oduced by the DI eve nt (hereafter, DI ejecta plume). 
Subaru Telescope observation (STO) by iKadono et al.l (120071 ) showed that the outer part of 
the DI ejecta plume (high-velocity ejecta) have a higher color temperature and lower infrared 
silicate-feature strengths than the inner part of the DI ejecta plume (low- velocity ejecta). 
They attributed the high-velocity ejecta to submicron-sized carbonaceous grains and the 
low-velocity ejecta to submicron-sized silicate grains. As a result, they proposed that the 
near-surface region of Tl is covered by a pristine crust of submicron-sized carbonaceous 
grains (hereafter, pristine carbon crust), which were created by cosmic rays in the trans- 
Ne ptunian region . In a ddition, they interpreted the existence of many craters observed 



by iThomas et al.l (120071 ) as evidence of a long surface retention age of the pristine carbon 
crust. This means that the nucleus of a comet does not form a dust mantle during its close 
approaches to the Sun. They also claim that the pristine carbon crust is very thin (several 
tens of cm) based on the ejection velocity and their estimate of the amount of carbonaceous 
grains. Finally, they concluded that even at a depth of only several tens of cm below the 
surface comet r naterials keep t heir p rimordial composition together with submicron-sized 



silicate grains. iFurusho et al.l (120071 ) supported this interpretation by measuring slightly 



higher linear polari zations for the high-velocity ejecta compared to the low-velocity ejecta. 
Harker et al.l (120071 ) also supported these interpretations from their infrared observations. 
These conclusions (hereafter, KFH model) completely contradict our conventional images of 
comet nuclei based on the standard model. 

Note that the difference in the color temperature between high- and low-velocity ejecta 
may simply reflect the amount of volat i le ma terials, because color temperatures depend 
on the volatile content (IMukai fc Mukailll984l ). It is also worth noting that the silicate- 
feature strength is not a good in dicator of an amount of silicate grains nor their sizes 
( lOkamoto. Mukai. fc Kozasal Il994j ) ; Large fluffy aggregates of tens of micrometers or larger 
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show strong silicate-features, while large compact aggrega tes of the same composition and 
mass show weaker silicate-features (IKolokolova et al.ll2007l ). Furthermore, the degree of lin- 
ear polarization for aggregates may slightly vary with their porosity and/or volatile content 
faimura et aUbooeh . jj^ 



In this study, we examine how well the standard model can explain the available obser- 
vational data of the DI event by calculating color temperature, silicate-feature strength, and 
polarization for aggregates in the framework of the standard model. We then discuss the ad- 
vantage of the stan dard model over the KFH model even when comparing wit h the results of 
Sugita et all (2005 ) who a nalyz ed the same STO d ata as Kadono et al. ( 2007 ) (hereafter, we 
refer iKadono et all (120071 ) and ISugita et al.l (120051 ) as to STO-K and STO-S, respectively). 
We also discuss the groundless interpretations of the long retention age of the Tl primor- 
dial surface and the thickness of its dust mantle based on comprehensive observational and 
theoretical considerations. 



AGGREGATE MODEL 



We calculate the color temperature T^, the silicate-feature strength Ssu and the degree 
of linear polarization Pi for the standard model as follows. In the standard model, large 
compact aggregates exist in a d ust mantle and large fl uffy aggregates are embedded with 
volatiles below the dust mantle (IKolokolova et al.l 120071 ). To describe lar ge aggregates, we 



use f ractal clusters of 2^^ identical monomers with a radius of 0.1 /im (IKolokolova et al. 



20071 ). Each monomer is assumed to have a concentric structure of an organic refractor y 
outer layer, a forsteritic inner layer, and an amorphous silicate core (lYamamoto et al.ll2007l ). 



The compact aggregates are assumed to have a fractal dimension D = 2.5, which 
corresponds to the value when the maximum compression of fluffy aggregates is attained 
(IWada et al.l 120071 ). Our calculation using the superposition T-matrix method (TMM) 
showed that the geometric alb edo of the compact aggregates is ~ 0.04 at visible wave- 
lengths (see iKimura et al.ll2003l ). and this value is identical to the albedo derived from an 
observation of the Tl surface (A'Hearn et al. 2005). We also calculate (3 ~ 0.3 (the ratio of 
solar radiation pressure to solar gravity) for the compact aggregates, using Mie theory with 



the o ptical constants deduced from the Maxwell-Garnett mixing rule (MG) (see lMukai et al. 



1992h: This is also consistent with f3 for high- velocity ejecta within the DI ejecta plume (e.g., 
Meech et al.ll2005h . 



Similarly, let large fluffy aggregates embedded wit h volatile materia ls below the dust 
mantle be presented by fractal clusters with D = 1.9 (iMukai et al.lll992l ). ISunshine et al. 



( 120071 ) found the 3 /im feature of solid H2O ice in their observations of the DI ejecta plume. 
Sublii nation of volatiles wi thin the ejecta seems to have lasted over the course of a few 
days ( Schleicher et aljbood ). We thus consider ongoing sublimation of volatiles in the fluffy 
aggregates and assume that a few percent of volatiles in volume still remain on the surface 
of the aggregates. (For a comparison, we also consider volatile-free fluffy aggregates). 

We then derive the absorption cross sections C\ for the compact a nd fluffy aggregates , 
using Mie theory with the optical constants deduced from the MG (e.g.. lMukai et al.lll992l ). 
The use of the MG for aggregates has b een examined based on the discrete dipole approx- 
imation (DDA) (IKolokolova et al.l 120071 ). From values of C\ at wavelengths A of 8.8, 10.5, 
and 12.4 /xm where STO-K observed, we calculate the color temperature as Tc = Cs.s/Ci2a 
and the sihcate-feature strength as = 2Cio.5/(C8.8 + Ci2a) — 1- We also calculate Pi at 
A = 0.8 /im and at a phase angle of the DI event for these aggregates. Since Pi does not 
strongly depend on the number of monomers, we perform numerical calculations with clus- 
ters of 256 monomers by using the DDA and the TMM. All the results are listed in Table 
1. 



3. STANDARD INTERPRETATION 



We first interpret the observational data by STO-K. Table 1 shows that the compact 
aggregate has higher Tc and lower Ssi than the volatile-coated fluffy aggregate. Impact crater- 
ing mechanism indicates that high- velocity ejecta are thrown from near-surface regi ons and 
low-v elocity ejecta contain primarily materials excavated from deeper regions (e.g.. iMelosh 
19891 ). Thus, a natural expectation in the standard model is that high- velocity ejecta in the 
DI ejecta plume would have higher Tc and lower Ss\ than the low-velocity ejecta. This agrees 
with the observational data by STO-K that show the high-velocity (~ 180 m s~^) ejecta to 
have higher Tc and lower S^i than the low- velocity ejecta (~ 100 m s~^). 



We next compare the polarimetric observation by lFurusho et al.l (120071 ) with Pi in Table 
1. Table 1 shows that the both aggregates have a positive polarization and that Pi for the 
compact aggregates (the high-velocity ejecta) is larger than that f or the fluffy agg r egate s (the 
low-velocity ejecta). This is consistent with the observation by iFurusho et al.l (120071 ) who 
showed that the high-velocity ejecta have a slightly higher Pi compared to the low-velocity 
ejecta. 

The STO data at time t = 3.5 hrs after the DI event are shown by STO-S as infrared 
spectra at various radial distances. We also interpret the STO-S data by estimating the radial 

-^8.8/-^12.4; 'S'gj = 2/io.5/(/8.8 + -^12.4) ~ 1; 



distributions of T^' and S'^j from the data, where T^,' 



I\ being the STO intensity at a wavelength A of 8.8, 10.5, or 12.4 yum, respectively. Figured] 
shows that both and S'^j at t = 3.5 hrs have the same peak position at ~ 3'.'0. We consider 
that the peaks represent low-velocity, volatile-free fluffy aggregates, because (1) the peak 
position is within the expected position (> 2'.'5) of the low-velocity ejecta at t = 3.5 hrs; (2) 
the peak position of the high- velocity ejecta at t = 2.0 hrs is already 2'.'5 (STO-K), so that 
the expected position aXt = 3.5 hrs is beyond ~ 4'.'0; (3) the volatiles in the fluffy aggregates 
would sublime gradually through a few hours, resulting in higher T^, (see Table 1). Thus, 
the STO data at t = 3.5 hrs is well explained by the standard model. 

In summary, the optical properties of the compact and fluffy aggregate s based on the 



standa rd model well account for the observational data used by STO-K, iFurusho et al. 



( 120071 ) ■ and STO-S. 



4. DISCUSSION 

In this section, we refute the following interpretations in the KFH model: (1) the 
presumption of the pristine carbon crust; (2) the long retention age of the pristine carbon 
crust; (3) the existence of cometary primordial materials at a depth of several tens of cm 
from the surface. 



4.1. Peculiar Pristine Carbon Crust 



In the KFH model, submicron-sized (0.3 /im radius) carbonaceous grains are responsible 
to the T'^ peak within the high-velocity ejecta (~ 180 m s~^). If this were correct, the 
carbonaceous grains in the high-velocity ejecta must have been decelerated significantly 
through t = 2-3.5 hrs, because the peak position of is ~ 2'.'5 at t = 2 hrs (STO-K) and 
proceeds to only ~ 3'.'0 at t = 3.5 hrs (Fig. [1]), corresponding to an average velocity of 60 m 
s~^. Furthermore, STO-K claimed that the high- velocity ejecta (> 160 m s~^) were observed 
at t = 6-9 hrs by iFurusho et al.l (120071 ) as highly polarized regions at > 6". Therefore, the 
KFH model suggests that the high- velocity (~ 180 m s~^) carbonaceous grains were suddenly 
decelerated to ~ 60 m s~^ during t = 2-3.5 hrs and then accelerated again to > 160 m s~^ 
after t = 3.5 hrs. We believe that no one could support such a dynamical behavior. A 
most likely explanation of the STO data is that the high- velocity ejecta expa nd with a 
nearly constant velocity over the first 20 hrs as observed by iMeech et al.l ( 120051 ) and that 
low-velocity ejecta is responsible to the peak att = 3.5 hrs. In addition, the observational 
constraint of /9 ~ 0.3 in the high-velocity ejecta completely rules out the predominance of 



0.3 /iin carbonaceous grains, which have (3 > 1 (IBurns et aLlll979l ). Therefore, the pristine 
carbon crust is a false interpretation, unless the above contradictions are resolved. 



4.2. Crater Formation on the Tl Surface 



STO-K interpret the existence of craters with radii up to ~ 200 m (IThomas et al.ll2007l ) 
as evidence for a long surface retention age of the pristine carbon crust, because they assumed 
that obs e rvable craters are produced when Tl was in the trans-Neptunian region. However, 
Kadond fll999l ) himself estimated that a comet collides Near-Earth Asteroids (NEAs) of 
Im every ~ 100 years during its approaches the Sun, forming a ~ 100 m sized crater. 
Short-period comets are also bombard ed in the main asteroid belt with higher collisional 
probabilities than those of NEAs (e.g., iHawkes fc Eatonll2004l ). The model calculation by 
Gronkowskil (120071 ) shows that the number of observable craters on Tl is around 10 for the 
steady state of collisional environment in the main asteroid belt and up to ^ 120 . Even if 



we use a different pa r amete r set with the asteroidal collision model by iGil-HuttonI (120001 ) 
and lHawkes fc EatonI (120041 ). the number of ~ 200m-sized craters produced at each orbit is 
~ 20. Thus, many craters with radii < 200 m can be formed during every approaches to th e 
Sun. Furthermore, sublimation processes increase crater sizes with time (IBritt et al.ll2004l ). 
suggesting that the craters observed on Tl may have been produced by smaller impactors 
with more frequent impact rates. We thus conclude that the existence of craters on Tl 
cannot be attributed to a long surface retention age. 



4.3. Thickness of the Tl Dust Mantle 



The sta ndard model predicts th e thickness of the dust mantle on Tl to be several meters 
or thinner (jPe Sanctis et al.l 120071 ). Following STO-K, we could estimate the maximum 
depth of the reg ion from which the compact aggregates are excavated, by using the Z-model 
(IMaxwelll 119771 ). This model allows us to estimate the maximum depth of an excavation 
flow (maximum excavation depth) for each ejecta with ejection velocity Ve- As a result, we 
might arrive at the dust mantle thickness of severa l tens of cm as STO-K, bu t we find several 
problems in their assumptions. According to lCroftI (Il980l ) and lMeloshl ( 1l989l ). the center of an 
excavation flow (flow center) in the Z-model is at some depth below the surface, which may 
be comparable to a projectile diameter even for oblique impacts. However, STO-K assume 
that the flow center is at the surface. This is not correct: the actual excavation depth is 
the sum of the depth of flow center (the projectile diameter, ~ 1 m) and the maximum 
excavation depth. In addition, an acceleration of the ejecta due to the gas plume observed 
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by STO-S was not taken into account. Then, the value of used by STO-K is overestimated 
and the resulting thickness is underestimated. Although STO-K claim that the thickness of 
several tens of cm is similar to that derived from the total mass of carbonaceous grains in 
the DI ejecta plume, the presumption of the pristine carbon crust is a false interpretation, 
invalidating their estimate of the total mass of the carbonaceous grains. As a conservative 
estimate, the DI results suggest that the thickness of the dust mantle is 1-2 m. 

This study was supported in part by T. Yoda, a Grant-in-Aid for Scientific Research 
from Japan Society for the promotion of Science, and a Grant-in-Aid for Scientific Research 
on Priority Areas from Ministry of Education, Culture, Sports, Science and Technology, 
Japan. 
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Fig. 1. — The radial distributions of the color temperature T^' and the silicate-feature S!,^ 
derived from the data of Sugita et al. (2005; STO-S) at t = 3.5 hrs. The radial distance is 
the value from the nucleus whose actual position is located somewhere in the point spreading 
function (~ 0'.'4) of STO (see, STO-S). 
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Table 1: Color temperature T^, silicate-feature strength Sgi, and linear polarization P\. 





T 

^ c 


'S'si 


Pi 


Compact aggregate 


1.94 


3.08 


0.211 


Volatile-coated fluffy aggregate 


1.62 


3.17 


0.167 


Volatile-free fluffy aggregate 


1.99 


3.51 


0.187 



